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ANALYTICAL DETERMINATION OF SPACE STATION RESPONSE TO CREW MOTION AND 
DESIGN OF SUSPENSION SYSTEM FOR MICROGRAVITY EXPERIMENTS 


by 

Frank C . Liu 

Professor of Mechanical Engineering 
The University of Alabama in Huntsville 
Huntsville, Alabama 


ABSTRACT 


The objective of this investigation is to make analytical deter- 
mination of the acceleration produced by crew motion in an orbiting 
space station and define design parameters for the suspension system of 
microgravity experiments. A simple structural model for simulation of 
the IOC space station is proposed. Mathematical formulation of this 
model provides the engineers a simple and direct tool for designing an 
effective suspension system. 
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1 . Introduction 


Some microgravity experiments to be performed on the IOC space 
station require an environment in which the acceleration level must be 
below 10"“ 5 g. Among the various sources of disturbances, crew motion is 
the severe one which can produce acceleration to a magnitude of 10~ 4 g. 
The objectives of this investigation are to define design characteristics 
for the suspension system and to find effective means for isolation of 
the experiment packages from disturbance. To achieve these goals, the 
following have been accomplished. 

1 • Analytical formulation of acceleration response of IOC to 
crew motion 


Assumed modes method and modal transformation are used to obtain 
mathematical solution for the normal coordinates to an input function 
due to crew motion. Acceleration response of any point in the space 
station can be analytically formulated. 

2. Structural modeling for IOC space station 

The finite element model for IOC space station, shown in Figure 
1, used by JFC [1] for vibration analysis generates a wide range and 
closely spaced spectrum of vibration frequencies (see Table 4.3. 3.3-3, 
reference 1). It makes it extremely difficult for mechanical design 
engineers to identify a particular mode which has dominant effect on the 
experiment package. A simple structural model which is made of a rigid 
main body (extension of keel frame which supports all the massive 
modules), a cantilever beam (the keel frame), and a transverse canti- 
lever beam (the solar boom) is proposed for simulation of the IOC. A 
complete analytical formulation for this model is made. This formula- 
tion provides the design engineers a simple and direct tool for deter- 
mination of acceleration and disturbing frequencies acting on the 
experiment package. Computations can be carried out by using a pocket 
calculator. This model can be easily improved by increasing the 
degrees-of-freedom. A simple computer program will do this work. 

3. Derivation of design characteristics for suspension system 

Consider that a microgravity experiment package is supported by 
a spring and a mass-spring-damper vibration absorber is attached to it. 
The acceleration magnification factor and the ratio of acceleration 
response of package to acceleration input of the support is then derived. 
The plots of this factor versus frequency ratio (IOC to mass-spring) 
provide engineers design parameters of the suspension system. 

4* Significant findings for effective isolation of disturbance 

The microgravity laboratory is 23.3/83.2 feet from the center 
of mass of IOC with/without the orbiter attached. It is shown by 



numerical study that the acceleration disturbance is greatly reduced when 
the orbiter is not present (i.e., microgravity lab is actually 59.9 feet 
farther away from the C.M. of IOC). This shows that the most favorable 
location beside the C.M. of IOC is the node of the fundamental elastic 
vibration mode of the IOC. 


2. Analytical Formulation 

2.1 Modal analysis of a structure system 

A brief presentation of modal analysis of a structure is given here. 
Let P be the coordinate of a generic mass point P in a structure and 
u(P,t) be its displacement which is expressed in the form [2] 

N 

u(P,t) “ £ *i^ q i^ ^ 2 ' 1 * 

i=l 

where i|>^(P) is an admissible function and q^(t) is the ith generalized 

coordinate of a set N. The equation of motion of the structure system 
in matrix form is 

[m] {q> + [k] {q> - {Q> (2.2) 

where [m] , [k] , and {q} are the generalized mass* stiffness* and force 
matrices, respectively. 

The kinetic and bending strain energies of the system are formu- 
lated from the following integrals: 

T - h f mptu(P.t)] 2 dP = m p '^ i 2 (P)dPq i 2 (2.3) 

B B 

V - kf EI[u” (P»t)] 2 dP - hEJ mpU.-tf)] 2 dP qi 2 (2.4) 

B B 

where mp is the mass density at P and El is the bending stiffness of the 

structure at P, the prime denotes partial differentiation with respect to 
spatial coordinates and the symbol B with the integral sign means inte- 
gration over the entire body B. From the above integrals* the elements 
of the [m] and [k] matrices are obtained respectively. 


m 


ij 


2 

3 Z T 

aq^qj 


j mpt|» i (P)<l>j 
B 


(P)dP 


(2.5) 
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EI *i (p >V ,(p > dp 


( 2 . 6 ) 


Denoting the force acting at P in the direction of u(P,t) by f (t), then 
the virtual work is p 


= Jf p (t)6udP f p <(' i (P)6q i dP 


(2.7) 


It follows that 


Q i " “ J Vi (p),lp 


( 2 . 8 ) 


Next, the natural frequencies of the system are determined from the 
determinant 


| [k] - ai 2 [m] | « 0 (2.9) 

and the eigenvectors (or modal columns) {<J>.} are the solution of the 
matrix equations, 1 

([k] - aj^tm]) {<1^} = {0} i = 1,2 N (2.10) 

The normal coordinates {ri} and the generalized coordinates are related 
by the transformation 

{q} = [#] {n} (2, 11) 

where the modal matrix 

[*] = [{+!>. U 2 h {* N }] . 

Applying Eq. (2.11) to Eq, (2.2) and making use of the property of 
generalized orthogonality of the eigenvector with respect to [m] and [k] , 
the equation of motion in normal coordinates is 

+ 55 N j (t)/M jj J = 1.2,. ...N (2.12) 

where 


M.. = {<(».} 


(2.12a) 



(2.12b) 


“j 2 “ { V T[kH V /M jj 

Nj (t) « {4 »j> T {Q> (2.12c) 

2.2 Structure response to disturbing force 

Consider that a disturbing force f (t) is applied to the structure 

a 

at P , from Eqs. (2.8) and (2.12c) one obtains 
a 

Y C > ■ E Vk <F a )f a (t > <2 ' U) 

k=l 

where is the kth element of the jth eigenvector. The response at 
point P e in the structure, u(P e ,t), can be obtained by using Eqs. (2.1), 
(2.11), and (2.13). It results in 


u(P e .t> - £ 

w>i (t > - E E *i ( y*ijV t) 


i 

i 

j 


-E 

|Y E Vi <p e>) 

(E ♦kjw)] 

(2.14) 

j 

L\ i / 

V k /J 



where (t) denotes the solution of the differential equation. 


*ti. +u». n. = f (t) with n.(0) = n. (0) = 0 
3 J J a J J 


(2.15) 


2.3 Model of crew motion and response function 

The motion of an astronaut inside a space module is started by push- 
ing one wall and motion is stopped by pushing the opposite wall. A 
simple mathematical model is suggested [1] as shown in Figure 2. 
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0 


(2.16) 


0 < t *- t. 


f a<‘> 


0 


C 1 < t < t 13 

C 13 < — t M 

fc 14 < 


The magnitude of f Q is 25 lbs and t^, t^ 3 » and tj^ are 1, 13, and 14 

seconds, respectively. These notations for time will be kept through 
the formulation for the sake that one may wish to change their magnitudes 
and secondly, that the equation involved will have appropriate units 
(use sina)t 13 rather than sinl3w where cu is in rad/s). The solution of 

equation (2.15) is readily obtained. Denoting that 

n(t) - (f 0 /u 3 t 1 )u(t) (2.17) 

the dimensionless displacement function U(t) is given by 


o 

A 

ft 

|A 

rt 

U(t) = cjt - sinmt 

(2.17a) 

' t < t 13 

U(t) * - sinwt + sinu(t - t^) + u)t^coso»(t - t^) 

(2.17b) 

c 13 * t - t 14 

U(t) = u»(t - tj^) - sinait + sinb)(t - tj) 

- sinui (t - t 13 ) + u)t^[cosh>(t - tj) 



+ cosw(t - t 13 )] 

(2.17c) 

rt 

♦— 

A 

rt 

U(t) = - sinajt + sina>(t - t^) - sinu(t - t^ 3 ) 
+ sinw(t - tj^) + a)tj[cosu)(t - t^) 



+ cosu(t - t ^ 3 ) ] 

(2. 17d) 


The function U(t) is continuous and has continuous first derivative 
(velocity); its second derivative is only piecewise continuous. Denoting 

n(t) = (f o /ut 1 )A(t) (2.18) 

and differentiating Eq. (2.17) twice, one obtains the dimensionless 
acceleration A(t), 
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0 < t <. A(t) = sinwt (2.19a) 

< t < t 13 A(t) * sintut - sinu)(t - tj) - utjcosu^t - tj) (2.19b) 

< t £ tj^ A(t) = sinwt - sinoi(t - t^) + sina)(t - t 13 ) 

- [cosw(t - tj) + cosw(t - tj 3 )J (2.19c) 

< t A(t) = sinuit - sini»)(t - t^) + sinu)(t - t^ 3 ) 

- sinu(t - t ^) - utj^ [cosw(t - t^ 

+ cosu(t - tj^)] (2.19d) 


The functions of U(t) and A(t) are plotted versus time for frequencies 
ranging from 0.1 to 0.4 as shown in Figures 3 and 4 respectively. Now, 
the displacement and acceleration response of point P g in the space 
station can be expressed in the form 


“<v*> - LIE 

♦uW) (E 

j ' i 


*<%.*> - e(e 

♦ij'Mv) (E 

i Vi 

' ' k 

The subscript "j" with 

U and A denotes 

ponding to u> = ok • 



(2.20a) 




(2.20b) 


2.4 Modeling of IOC space station 

The structure of the IOC space station may be treated as a struc- 
tural system having three elements. The main body is a frame structure 
which supports all the massive members, the vertical and horizontal HAB 
modules, the vertical and horizontal LAB modules, logistic and common 
modules, and above all, the orb iter. Attached to the main body is a 
keel frame structure 296 feet in length which supports an antenna system 
at its other end and a transverse boom at a distance 165.5 feet from the 
main body. The third member of the system is a transverse boom which 
is a frame structure 264 feet long. Its main purpose is to carry eight 
solar arrays and power system radiators. 


» 
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A simple structural model for the IOC proposed consists of a rigid 
main body which supports a cantilever beam (keel frame) and a transverse 
cantilever beam (solar boom) mounted on the keel. As shown in Table 1 
in Appendix A, the modules, equipment, solar arrays, fuel tanks, orbiter, 
etc., are treated as concentrated masses. However, the rotational 
moments of inertia of the orbiter and solar arrays must be included in 
forming the mass matrix. 

• 

Motions of the IOC in X— Z and Y— Z planes will be treated separately. 
In each plane the proposed IOC has four degrees-of-freedom, namely, 
rigid-body translation, rigid-body rotation, and one bending mode for 
each cantilever beam. Thus, the corresponding admissible functions are 


CP) = 1 (2.21a) 

= (2.21b) 

4- 3 (P) = z 2 - z 3 /3 (2.21c) 

^ 4 (P) = y 2 ~ y 3 /3 (2.21d) 


where z * z/£^ and y = y/£ g in which and i, g are the length of the keel 
(296 ft) and solar boom (132 ft), respectively. 


2.5 Modal analysis of IOC model 

Based on Eq. (2.21) the mass and stiffness matrices for the model 
are formulated in Appendices A and B. The matrix equation of free 
motion of the IOC is 



( 2 . 22 ) 


One may eliminate the rigid-body translation q^ from the system by solv- 
ing q^ in terms of the rest of q's from the first equation of Eq. (2.22) 

and substituting it into the remaining equations. A further simplification 
can be made by disregarding the small coupling effect of the rigid-body 
rotation q^ and the elastic modes q^ and q^ due to gravity gradient 

torque. Thus, one may eliminate both q« and q« from the system and obtain 
the following: 1 1 
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(2.23) 



(2.24) 


where 


[T] = - 

m ll 

m 12 

-1 

m 13 

m 14 1 


_ m 21 

m 22_ 


m 

CM 

a 

i 

m 24 J 



- 


- 



r i 

m 33 

“34 

s 

m 33 

“34 

+ 

“31 

“32 

“4 3 

-F> 

1 


“43 

“44. 


.“41 

CM 

e* 


The two elastic frequencies of the IOC can be written out directly 
in the form 


+ Xm^ 


0 ), 


m 33 k 44 + m 44 k 33 + /(m 33 k 44 m 44 k 33 ) + 4m 34 m 43 k 33 k 44 


1,2 

and the eigenvector is 


2(m 33 m^ - m 34 


(2.25) 



T ll (m 34 

“i 2 > 

+ 

T 12 

( k 33 

- m 33 

uk 2) 


T 21 ( “34 

u i 2) 

+ 

T 22 

( k 33 

' “33 

u, i 2) 

{♦ ± } = 




2 






m, 

33 

“i 





_ 

k 33 

- 

“33 

2 

u i 




(2.26) 
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2.6 Rigid-body librational motion of IOC 

The equations of librational motion of an orbiting body which has 
its principal axes parallel to the orbital axes subjected to disturbing 
torque M can be written directly in the form [3] 


l A + 3 “oS - V\ - « x 

Vy + 3 “o 2<I x - ■ «y 


(2.27) 

(2.28) 


where the I's are the moments of inertia about their respective principal 
axes through the center of mass. The orbital frequency of a circular 
orbit is 


u 

o 



rad./s 


(2.29) 


where R is the orbital radius and the gravitational constant 

y - 1.407 x 10 16 ft 3 /s 2 
= 3.986 x 10 14 m 3 /s 2 

Thus, Eqs. (2.27) and (2.28) yield the librational frequencies 

“x = “o /3(I y " V /r x (2.30a) 

w y = % /3 < x x " (2.30b) 

2.7 Response of IOC to crew motion due to rigid-body modes 

As shown in Figure 6, due to the arrangement of the modules, crew 
motion will create a disturbing force either in x-direction (motion in 
horizontal modules) or z-direction (motion in vertical modules). Hence, 
the torque M x is negligible in comparison with M which has the magni- 
tudes, y 


M 

y 


f z 

x a 


,f X 
z a 


(motion in horizontal modules) 
(motion in vertical modules) 


(2.31) 


As indicated by Eq. (2.30), the librational frequency is approximately 
equal to JT times the orbital frequency (I x a I , « I x » see Table 1). 

The crew kicking motion is completed in a time interval of 1 second which 
is very short in comparison with the period of librational motion of 3,300 
second. This means that the torque can be treated as an impulse torque. 
The action of an impulse torque is equivalent to give the space station 
an initial angular velocity, i.e., 

t 

My = f My(t)dt - e y (0) Iy (2.32) 

o 


The solution of Eq. (2.28b) is simply 


9 » (M /u) )sinw t (2.33) 

y 77 y 

Using the mathematical model for crew motion given by Eq. (2.16), one 
obtains the response at P ft 

“(P e ,t) - Z e 0 y 

Z Z f t,/« I )sinu> t (motion in horizontal module) 
e a o 1 y y' .y 

Z X f t,/w I )sinu t (motion in vertical module) 
e a o 1 y y y 

(2.34) 


The magnitude of acceleration produced by crew motion at P^ is 


[ Z Z f hi t,/2I (motion in horizontal module) 

e a o y 1 y (2 . 35) 

Z X f (D t,/2I (motion in vertical module) 
e a o y 1 y 



3. Numerical Results of Acceleration Due to Crew Motion 

Using data given in Table 4. 3. 3. 4-2 to 4. 3. 3. 4-4 of Reference 1, 
Table 3 is formed for the formulation of the mass matrix for IOC space 
station. Note that the total weight given by Tables 4. 3. 3. 4-2 and 
4. 3. 3. 4-4 [1] is 77,600 lbs heavier than that given by Table 2. In an 
effort to match the total weight given in Table 2, some of the weights 
are not included. There is a significant difference on the location of 
the center of mass between the present model to that given in Table 2, 
as shown in the following: 
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Table 2 


Present Model 


Difference 


With Orbiter 68.3 ft 113.8 ft 14.4 ft 

Without Orbiter 128.2 59.5 8.8 

Consider that the distance from the reference point to the micro- 
gravity experiment is equal that of the location of crew motion, i.e., 
z e = z c = ”55 ft. The magnitude of acceleration on the microgravity 

experiment due to crew motion as given by Eq. (2.20b) is 

a e = t*lj + VVV’ 2 [A j (t)I ma* 3 ‘ U2 <3a) 

The value of [A^ (t) ] may can be estimated from curves given in Figure 4 

for a given value of w. The summation is omitted so that a is calcu- 

e 

lated for each mode. Based on the librational motion approach, one has 
from Eq. (2.35) 

2 

a * Z f ai t. /2I (2 2 ) 

e e o y 1 y w •*-) 

The numerical resulted obtained are summarized in Table 4. 

It is important to note the following: 

(1) No direct comparison can be made on the frequencies obtained 
to that given in Table 4.3.3. 3-3 [1] due to the difference of inertia 
properties of the models, and furthermore, the JSC model has no distinct 
fundamental mode that can be singled out. 

(2) The acceleration level obtained here is about one order smaller 
than that given by Table 4.3.3.5-15 [1]. This is due to the fact that 
the inertia data given by Table 4.3. 3. 5-6 [1] is about 1/3 of that given 
by Table 4. 3. 3. 5-6 [1] (without orbiter). In addition, the disturbance 
torque given by Table 4.3. 3.5-7 [1] is more than 2 times the value used 
here. Thus, the magnitude of accelerations presented in Table 5 are 
reasonable. 

(3) Motion of the space station in Y— Z plane will occur if the 
disturbing force is in the direction parallel to the solar boom. 

(4) The acceleration given by the librational motion is 10“^ of 
that given by elastic motion. This, due to the frequency of librational 

-3 

motion, is only 10 of the frequency of elastic motion. 


4. Design Characteristics of Suspension System 


4.1 Acceleration response of a vibration absorber system 

Consider that a microgravity experiment package is mounted on the 
laboratory module structure which has motion u Q (t) as a result of crew 

motion or other disturbance. As shown in Figure 3, u Q (t) is approximately 

a harmonic. It is required to design a suspension system which can 
effectively reduce this disturbance over some frequency range. Vibration 
can be effectively reduced by using a vibration absorber [4] which is a 
mass-spring— damper system attached to the main mass-spring system as 
shown in Figure 7. 

Denoting u q , u^, u 2 » the absolute displacement of the structure, 

main mass, and absorber mass, respectively, the equations of motion of 
the system are 

In l'*l + c “l + ^1 + ^2^ u l ~ c ^2 ~ ^2 U 2 * kjU Q (t) C4.1a) 

m 2^2 + c “2 + ^2 U 2 ” C1 *l ~ ^2 U 1 * ® 


Letting the input u Q (t) be a harmonic disturbance, one may put 


U Q (t) = U o e iuit , Uj - Uje 1 ^ , and u 2 * U 2 e iU,L ' i - /-T (4.2) 


iut 


iut 


where U 1 and U 2 are complex quantities that can be determined from the 
matric equation 


_ o 

(k^ + k 2 - m^u + icu) 

- (k 2 + icu) 


1 

H 

:=> 

i 


k.U " 
1 0 

- (k 2 + icu) 

o 

(k 2 - m 2 u + icu) 


L U 2- 


0 


(4.3) 


If one wishes to determine acceleration response rather than displace- 
ment , set 

Uj * A^e^ u)t with A^ * - u^B^ (4.4) 

Now, let M denote the magnification factor of acceleration, the ratio of 

3 . 

A^ to input acceleration. 
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(4.5) 


« a * UjUV 2 

^ rst » introduce the following dimensionless parameters: 

k = kj/k 2 M = m 2 /m 1 x, = c /2 /m 2 k 2 r = to/to^ ( 4 . 6 ) 

It can be shown that 

M a = { [(1 " H kr2 > 2 + 4pkc 2 r 2 ]/D | 1/2 ( 4 . 7 ) 

where 


D - |l - [1 + (l+k)y]r 2 + ykr 2 | 2 + 4yk£ 2 r 2 [1 + (l+y)r 2 ] 2 
4 . 2 Design considerations 

The magnitude of the magnification factor of acceleration depends on 
four parameters: 

k, the spring ratio (main spring /absorber spring) , 

U» the mass ratio (absorber mass/main mass), 

damping factor (damping coefficient/critical damping), 

r, frequency ratio (space station frequency /natural frequency of 
main mass-spring system) . 

To plot versus frequency ratio squared as shown in Figure 8 , six sets 

of curves are illustrated for ?= 0.2. The first three sets are for 
fixed values of y = 0.01, 0.025, and 0.05, respectively, with various 
values of k. The next three sets are for fixed values of k = 10, 15, 
and 20, respectively, with various values of y. All these curves have 

one common characteristic, which can be effectively reduced for r 2 > 2 . 

Since the frequency of the IOC space station in elastic vibration is about 
ra /s, it requires that k^ > 2W^/g. For a microgravity experiment 

package of 1,000 lbs, the spring constant of the suspension must be less 

than 5 lb/in. For the frequency range 0.5 < r 2 < 2, can be made less 

than 0.5 by a proper combination of y and k. 
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5. Conclusions 


A simple structural model for simulation of IOC space station has 
been presented and formulation of this simple model provides engineers 
a simple and direct method for computing the fundamental frequencies of 
the space station and determining the magnitude of acceleration at any 
point produced by crew motion. Acceleration response of a mass-spring- 
absorber system to a moving support is also formulated. Design engineers 
can use plots of acceleration magnification factor versus frequency 
ratio squared to determine design parameters for the suspension system 
of the microgravity experiment package. 

The following are some significant findings: 

(1) The acceleration due to librational rigid-body motion is 3 
orders smaller than that due to elastic bending motion. 

(2) The frequency of librational motion is 1.8 x 10 rad/s, which 
is approximately vT times the orbital frequency. 

(3) The frequencies of elastic motion of the simple model are in 
the range 0.6 to 1.5 rad/s. 

(4) Only the fundamental bending mode has dominant contribution to 
the acceleration, therefore, a simple model is adequate. 

(5) An effective suspension system can reduce the acceleration to 
1/4 of its magnitude, at most. 

(6) The ideal location of the microgravity lab does not have to be 
near the center of mass of the space station. 

(7) The most effective means to eliminate acceleration is to have 
the experiment module near the nodal point of the fundamental bending 
mode, as illustrated in Figure 9. This means that the factor (<J> 1X + 

d) z /a ) in Eq. (4.1) becomes very small. Example: The IOC without 

21 e k 

orbiter has moved the lab module 59.9 ft further away from the center of 
mass, but the acceleration is reduced 3 orders smaller than the IOC with 
orbiter. 

(8) It is favorable to perform the microgravity experiments when 
the orbiter is not present. 

(9) It is possible that by rearranging some massive elements, a 
minimum value of (<f>^ + ^21 Z e^k^^ll Can reac ^ ie< ^' However, this 

has to be done by trial and error method. 
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Notations: 


Appendix A. Formulation of Mass Matrix 


I , I 

sx’ sy 


^ox * ^ 0 y moment of inertia of orbiter about x- and y-axis through 
C.M. of orbiter, respectively 

I s x J ^y moment of inertia of solar arrays about x- and y-axis 
through solar boom at the attachment, respectively 

I UB moment of inertia of upper boom about x-axis through end 

of keel 

moment of inertia of lower boom about x-axis through 
attachment 

length of keel frame 
£ length of solar boom 

b 

mass of keel frame alone 

m, kth concentrated mass attached to keel at distance z, 

k 

mass of solar boom frame alone 

m s sth concentrated mass attached to solar boom at y g 

m r rth concentrated mass attached to rigid main body at z r 

m mass of orbiter 

u x » U y displacement in x and y direction respectively 

y s ~ y s /* s . coordinate of concentrated mass attached to solar boom 

\ ~ z i c ^ 2 'k* coordi nate of concentrated mass attached to keel 

z 0 coordinate of orbiter 

subscripts & 

su mm ation index k(keel) , r(rigid main body), s(solar boom) 


u , u 
x’ \ 


A-l. Motion of IOC in X-Z Plane 


The kinetic energy of rigid main body is 
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t r ■ 1/2 L m A 2 (z r ) 

r 

= 1/2 ^ m r^l^ z r^l ^2^ z r^2^ "*" ^oy^2^ z o^ ^2^ 

r 

The kinetic energy of keel structure is 

T K = 1/2 21 W'k^l + ^2 (z k ) ^2 + +3 (z k )4 3 lZ 
k 

+ 1/2 + * 2 ( 2)42 + ^ 3 C z ) 4 3 ] 2 dz /\ 

+ 1/2 Z I sy l ' l '2 (z s )4 2 + *3 (z s^3 l2 
s 

The kinetic energy of solar boom is 

Tg = 1/2 M g J[* 1 (z s )q 1 + * 2 (z s )q 2 + VV'U^s^s 

* 2 

+ 1/2 £ m s [i|; 1 (z s )q 1 + ^ 2 (z s )4 2 + * 4 (y s )4 4 l 

S 

+ 1/2 L ww -^ 2 

s 


(A-l) 


<A-2) 


(A- 3) 


By using Eq. (2.5) the elements of the mass matrix are obtained as follows: 

m ll ■ “total " M r +M k + £ m k + M s + 2I m s 

k s 

m 12 = Z m r z r + 1/2 “k + L Yk + (M s + L m s )z s 
r k s 

m 13 = 1/4 “k + 21 “k (i k 2 “ 1/3 \ 3) + (M s + 21 m s )(z s 2 ‘ 1/3 Z s 3) 
k s 
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1/4 M s + £ m s (y/ - 1/3 y/) 

s 

E Vr 2 + 1/3 \ + E Vk 2 + < H s + E °s )T s 2 + E r sy /l k 2 + V*P 

r k s s 

m\ + £ "A 3 ' 1 - 1/3 + (“s + E "s )(1 - »» V r s 3 

k s 

+ <E v v z - vV 

S 

Z YU, / 

s 14 

(? s 2 - 1/3 7 3 ) « u 

( 11 / 105 )^+^ “k ( \ 2 - 1/3 \ 3 > 2 + <«s + E m s )( V - 1/3 r s 3 > 2 

k s 

+ (Y I ) ( 2z - z 2 ) 2 /JL 2 
*** sy v s s k 

s 

(11/105)M S + £ m s (y s 2 - 1/3 y, 3 ) 2 + «„lly sl (2 - y s[ ) 2 

s 

+ [ y S 2 (2 - y S 2 )Z}n s 2 

Motion in Y-Z Plane 

The kinetic energy of rigid main body is 

T R * 1/2 £ VW’l + ^2 (z r^2^ 2 + 1/2 I ox t ' < '2 (z o ) ^2 ] 2 
r 

+ 1/2 I UB [^(z UB )^) 2 


(A-4) 



The kinetic energy of the keel structure is 


T K * L + *2* z k^2 + ’*’3^ z k^3^ 2 

k 

+ 1/2 «k /(♦,’! + <l» 2 ( z k^2 + *3^ z k^3^ 2dz k^k 


+ 1/2 I UB [^(z s )q 2 + *3( z s )q 2 r 


(A-5) 


The kinetic energy of the solar boom structure is 


1/2 M s J 1*^! 

s 


+ 'l'2 (z s ) ^2 + ^s^^s^s 


+ 1/2 m s j {[^'q 2 + ^(z s )q 3 ]y s - V^V** 

S 

+ 1/2 22 <® 8 t*l*l + ^ 2 ^ z s^2 + *3^ z s^3^ 2 


• 1 2^ 


+ t<j^q 2 + ^3(z s )q3]y s - > 


+ 1/2 Z 1 ** i *2*2 + *3 (z s )4 3 ■ 


(A-6) 


Applying Eq. (2.5) results in the mass matrix for motion in Y-Z plane: 


m u - Mr - (m u ) x2 


m 12 “ *®12*xz 


m 13 “ ^ m 13^xz 
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0 



"23 ‘ <“23 ) X2 + *2*3< Z 8 H £ <sx * V * Ev/ + l,3 'V. ! l 

8 S 

+ 1 W*2»j < ‘k > 

"24 * ‘[M V. + £ V.W * £ 1 «**<y,>]*2 

8 8 J 

"33 - <"33>*z + f£ V. 2 + £ «„ - V] [ *3 (z s” 2 + 1 W f *3< t k>l 2 

"34 ■ *[l5 Vs + £ "sVV*. + £ 

8 8 

m 44 •" (® 44 )*z 

The subscript "xz" denotes element of mass matrix of motion in X-Z plane. 

A-3. Mass and inertia properties for IOC model 

Table 3 is formed based on data given by Table 4. 3. 3. 4-2 to 4. 3. 3. 4-4 [1] 
for the purpose of formulation of the mass matrix. 
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Appendix B. Determination of Stiffness Matrix 


B-l. Bending strain energy 

The bending strain energy of keel and solar boom structure is given by 

A 2 A 2 

V - 1/2 J (EI) k [*5q 3 ] Z dz k + 1/2 J (EI) 8 [*J(y 8 )q 4 ] Z dy 8 

'o o 

Applying Eq. (2.6), one obtains the following non-zero elements of the 
stiffness matrix 


k 33 - 4(EI) k /3t k 3 - 4<EI),/3t s 3 


B-2. Moment due to gravity gradient 

Using the formula given by Reference 3, the following are obtained: 
Motion in X-Z Plane 

M y - - 3 (h/R 3 )(I x - I z )q 2 , k 22 » 3(y/R 3 )(I x - I z > 

in Y-Z Plane 

3(p/R 3 )(I y - I z )q 2 , k 22 - 3(p/R 3 ) (I y - I z ) 

where the moments of inertia are about the axes through the c.m. of the 
IOC space station. 


Motion 


M 
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Appendix C. Numerical Results 


C-l. Motion in X-Z plane (IOC with orbiter attached) 


18900 

-2758 

1096 

171.6 

<d 2 - 1.0431 

-2758 

3168 

831.9 

96.0 

1 

1096 

831.9 

483.7 

43.7 

uj 2 2 - 6.64 

171.6 

96.0 

43.7 

122.9 



1.1135 


.2770 

.9546 

U 2 ) - 

.1165 

-10.118 

-6.49 

.1639 


-411.87 


M u - 29157 
M 22 - 1488 


C-2. Motion in X-Z plane (IOC without orbiter attached) 



11590 

376.2 

1096 

171.6 

a) 2 - 2.0243 

376.2 

1829.6 

831.9 

96.0 

1 

1096 

831.9 

483.7 

43.7 

w 2 2 - 7.436 

171.6 

96.0 

43.7 

115.2 




.1936 

1.054 

-2.429 

.1095 


U 2 > 


.9485 

4.769 

-8.923 

-17.603 


M n - 191.5 
M 22 - 32169 


C-3. Motion in Y-Z plane (IOC with orbiter attached) 


[m] 


18900 

-2758 


-2758 

3298 


{*!> 


1096 

872.3 


0 

-92.5 


1096 

872.3 

513.5 

-74.5 


0 

-92.5 

-74.5 

122.9 


.3673 

1.194 

-3.332 

.1674 


U 2 ) ” 

.3332 

.9829 

-3.274 

-5.779 


u^ 2 - .8022 
w 2 2 - 7.881 

M u - 959.8 
M 22 - 3312 
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C-4. Motion in Y-Z plane (IOC without orbiter attached) 


11590 

369.8 

1096 

0 

369.8 

3298 

872.3 

-92.5 

1096 

872.3 

513.5 

-74.5 

0 

-92.5 

-74.5 

11 


oij 2 - .4301 
u> 2 2 - 7.046 


{*!> 


.1906 

.510 

-2.191 


-1.716 


{* 2 > 


.3221 

.5438 

-3.589 


-13.17 


M u - 927.5 
M 22 - 18582 


C-5 Computation of accelerations 

Using the data given above, Eq. (3.1), and with the aid of Figure 4 
for the value of [A(t) ] mav . » acceleration of the experiment package can 

be calculated. The results are shown in Table 4. 
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Table 1. Properties of Analytical Reference Configuration 
Space Station Model (Table 4. 3. 3. 3-1 [1]) 


Component 


Bending Torsional Mass 
Stiffness Stiffness Length 
(f t-lbs-ft) (Slugs/ft) 


Bending Torsional 
Strength Strength 
(ft-lbs) 


Booms and Keels 1.31 E+9 3.18 E+8 

30- Inch Astro 3.13 E+6 2.08 E+5 

Mast 


0.25 35,000 15,000 

2.3 3,480 208 


Table 

2. Inertial Properties of Analytical Model 



(Table 4. 3. 3. 3-2 [1]) 



Weight 

C. G. Coordinates 

Moments of Inertia 

Case 

(lbs) 

(ft) 

(lb-f t-sec^) 

Without 

269,000 

(1.1,0,84.2) 

(8.63 E+7, 

Payloads 



7.82 E+7, 

and Orb iter 



1.20 E+7) 

With 

373,200 

(-1.1,0,128.2) 

(2.06 E+8, 

Payloads 



1.98 E+8, 

Only 



1.45 E+7) 

With 

508,800 

(5.7,0,34.4) 

(1.45 E+8, 

Orbiter 



2.37 E+8, 

Only 



1.35 E+7) 

With 

608,600 

(3.59,0,68.3) 

(3.21 E+8, 

Payloads 



3.14 E+8, 

and Orbiter 



1.62 E+7) 



Table 3. Masses and Inertia Properties Used for IOC Model 


r 

Masses attached to main body 


z r ft 

z r 

ft 

W r lbs 

i 

Fuel tank and gases 


94.5 

- 5.5 

17,002 

2 

Module radiators 


82.5 

- 17.5 

3,000 

3 

Five-bay platform 


76.5 

- 23.5 

1,403 

4 

COM 1203 


64.0 

- 26.0 

11,000 

5 

Logistics module 


71.5 

- 28.5 

37,823 1 

6 

OMV and kits 


67.5 

- 32.5 

23,750 

7 

Horizontal lab 


45.0 

- 55.0 

54,295 

8 

Vertical lab 


34.5 

- 65.5 

27,067 

9 

Keel extension 


27.0 

- 73.0 

970 

10 

Vertical lab 


19.5 

- 80.5 

47,709 

11 

Lower boom 


13.5 

- 86.5 

728 

12 

Horizontal lab, SAA0207 and SAA0201 

9.0 

- 91.0 

37,089 

13 

Orb iter 


-26.66 

-126.66 

235,400 

s 

Masses attached to solar boom 

Z ft 

z ft 

V ft 

W lbs 



s 

s 

' s 


s I 


Solar boom structure 

265.5 

165.5 

0-132 

2,345 

1 

Power system radiators 

265.5 

165.5 

54.0 

750 

2 

TDM 2010 

265.5 

165.5 

63.4 

1,540 

3 

4 inboard solar arrays 

265.5 

165.5 

78.0 

4,787 

4 

4 outboard solar arrays 

265.5 

165.5 

132 

4,787 

• k 

Mass attached to keel 


Z k ft 

z k £t 

W fc lbs 


Upper and lower keel structure 


100-396 

0-296 

2,504 

i 

Remote manipulator 


162.0 

62. 

.0 

2,000 

2 

Refuel attachment, tanks and tools 

107.5 

7. 

.5 

4,625 

3 

TDM 2570 


210.5 

110. 

.5 

2,000 

4 

Instruments and storage shelter 


212.5 

112. 

.5 

4,625 

5 

Storage boxes and tools 


272.5 

172. 

,5 

9.850 

6 

Service attachments 


290.5 

190. 

.5 

3,750 

7 

TDM 2560 


295.5 

195. 

5 

7,055 

8 

Satellite 


324.5 

224. 

5 

20,000 

9 

Upper boom and antenna system 


396.0 

296. 

0 

17,734 


Rotational Moment of Inertia 


I ft-lb-s 2 

I ft-lb-s 2 


Orbiter (about c.m. of orbiter) 


7 x 

LO 6 

8 

x 10 6 


Solar array (about attachment) each 






Parallel to Nadir 


101,000 

107,300 


Normal 


2,790 

107,300 


Upper boom with antenna system 


158,800 

small 


Module radiators (about attachment) 

583,500 


4,470 


Power system radiators (about attachment) 

750 

25,120 


Note: Z is measured from the bottom end of the keel extension and z is 

measured from the joint of keel and keel extension 
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Table 4. Numerical Results for IOC Model 




Motion in 

X-Z Plane 

Motion in 

Y-Z Plane 


u) and 

With 

Without 

With 

Without 

Motion 

a 

e 

Orbiter 

Orbiter 

Orbiter 

Orbiter 


03 rad/s 

1.021 

1.423 

0.896 

0.656 


a 

1.14 g x 10“ 5 

1.44 g x 10 8 

0.86 g x 10' 5 

1.17 g x 10 -5 

Bending 

e 

u) rad/s 

2.577 

2.730 

2.810 

2.654 


a 

e 

2.26 g x 10" 6 

1.38 g x 10 ^ 

0.66 g x 10" 5 

2.9 g x 10" 6 

Libra. 

03 rad/s 

1.89 x 10~ 3 

1.89 x 10 -3 

1.85 x 10“ 3 

1.81 x 10 -3 

Motion 

a 

e 

0.254 g x 10~ 8 

5.13 g x 10 -8 

0.24 g x 10 -8 

4.72 x 10“ 8 



